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ABSTRACT 

A long standing problem in weak lensing is about how to construct cosmic shear 
estimators from galaxy images. Conventional methods average over a single quantity 
per galaxy to estimate each shear component. We show that any such shear estimators 
must reduce to a highly nonlinear form when the galaxy image is described by three 
parameters (pure ellipse), even in the absence of the point spread function (PSF). In 
the presence of the PSF, we argue that this class of shear estimators do not likely exist. 
Alternatively, we propose a new way of measuring the cosmic shear: instead of averag- 
ing over a single value from each galaxy, we average over two numbers, and then take 
the ratio to estimate the shear component. In particular, the two numbers correspond 
to the numerator and denominators which generate the quadrupole moments of the 
galaxy image in Fourier space, as proposed in Zhang (2008). This yields a statistically 
unbiased estimate of the shear component. Consequently, measurements of the n-point 
spatial correlations of the shear fields should also be modified: one needs to take the 
ratio of two correlation functions to get the desired, unbiased shear correlation. 
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1 INTRODUCTION 

Weak gravitational lensing refers to the weak and system- 
atic shape distortions of background source images (galax- 
ies, CMB, etc.) by the foreground inhomogeneous density 
distributions on cosmological scales. Since this effect only 
involves gravity, it has been widely used as a direct probe 
of matter density fluc tuations of our Universe (see, e.g., 
iHoekstra fc JainlboOSi for a recent review). 

The weak lensing effect can only be probed statistically 
due to the fact that the intrinsic projected shape of each 
galaxy is always somewhat random and anisotropic. A cen- 
tral theme in the study of weak lensing is to find unbiased 
cosmic shear estimators on galaxy images. This is indeed 
very challenging because the shape distortion due to weak 
lensing is generally much weaker than the intrinsic variations 
of the galaxy shapes. 

In the early stage of this field, most of the 
work focused on issues regarding the use of the 
quadrupole moments of a galaxy image as a shear 
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estimator 


Tvson et al.l Il990l: iBonnet & Mellied Il995l: 


Kaiser et al. 


1 19951: 


LuDuino & Kaiser 'l997|: Hoekstra et al.l 


19981: iRhodes et al 


12000.: Kaiser 20001. Ever since then, a 



number of other shear estimators have been considered in 
the literature, includi ng moments d efined by a certain set of 
orthogonal functions ([Bridle et a~! !2001: Bernstein Sz Jarvi3 
2OO2I: iRefreeier fc B acon' '20031: iMassev fc Refregieil l2005l : 
Nakaiima fc Bernstein. 2007. ), t he spatial derivativ es of the 



galaxy surface brightness field (|Zhanell200il2010al ). etc.. 

Conventionally, for each shear component, the shear es- 
timator is simply one number derived from a galaxy image, 
whose statistical mean is supposed to be equal to the true 
shear value, provided that the intrinsic galaxy image is sta- 
tistically isotropic. Unfortunately, even in the absence of the 
PSF, we show that such shear estimators at least do not ex- 
ist in simple forms, making them hard to use in practice for 
precise shear measurements {e.g., in the presence of noise). 
In the presence of the PSF, we argue that such shear es- 
timators do not likely exist. We give reasons for the above 
statements in iJJ] (The readers who are just interested in our 
new way of measuring shears may skip this section.) 

In ^ we present a new form of shear estimators: in- 
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stead of having only one number from each galaxy image 
for each shear component, one can keep two numbers, and 
use the ratio of their averages over many galaxies to accu- 
rately measure the cosmic shear. We find that this new way 
of measuring the shear can be easily implemented by using 
the method of Zhang (2008) (Z08 hereafter). The new type 
of shear estimator requires weak lensing statistics such as the 
n-point correlation functions of the shear field to be carried 
out in a slightly unusual way, but with little additional cost. 
This is discussed in ^ In ij5l we give numerical examples. 
Finally, we summarize in ^ 



in terms of the convergence k = ($11 + $22)72 and the two 
shear components 71 = ($11 — $22)/2 and 72 ~ ^12- 

The quadrupole moments of the lensed galaxy image 
are defined as follows: 



Qi. 



d XXiXjfLix) 



(2) 



where the origin of the coordinates has been chosen to be 
the center of the light, i.e., 



d^xxfLix) — 



(3) 



2 CONVENTIONAL SHEAR ESTIMATORS 

Conventional shear estimators are defined as a class of shear 
estimators, which average over a single quantity per galaxy 
to estimate each shear component. Most of the existing shear 
measurement methods belong to this class. For example, 
the method by Kaiser et al. (1995) and its extensions basi- 
cally use the quadrupole moments of each Gaussian-Profile- 
Weighted galaxy image t o measure the shear components; 
in the shapelets method (|Refregierll2003l ). the value of each 
shear component is estimated from best-fitting a shapelets 
model to each observed galaxy image; the method of Bern- 
stein and Jarvis (2002) evaluate the shear components by fit- 
ting each galaxy shape (also the PSF) with a se ries of orthog- 
onal 2D Gaussian-based functions (see, e.g.. iMassev et ahl 
I2OO7I for more examples). A common feature of these meth- 
ods is that they all generate one quantity per galaxy for each 
shear component. 

In i]2.1l we start our discussion with the well-known 
examples of shear estimators consisting of quadrupole mo- 
ments of galaxy images, and show what the issues are. In 
i]2.2l we show that, even in the absence of the PSF, any 
conventional shear estimator has to reduce to a highly non- 
linear form, making it hard to use in practice. In H2.31 we 
provide arguments as to why we think that conventional 
shear estimators do not likely exist when a PSF is present. 

2.1 A Review of the Problem 

The use of galaxy quadrupole moments as shear estimators 
has been a central topic in weak lensing for many years. It is 
therefore easier to start our discussion with the quadrupole 
moments. To present the issues clearly, let us first consider 
the case without the PSF or any photon noise. For conve- 
nience, we use (2:1, 2:2) or (x,y) instead of {6x,0y) for coor- 
dinates in 2D in this paper. 

Suppose that the surface brightness field of the lensed 
galaxy image is fhix^) on the image plane, and that of the 
original (pre-lensing) galaxy image is fs{x^) on the source 
plane, where x^ and x^ are the position angles on the im- 
age and source planes, respectively. We have the following 
relations: 



/i(x^) = /s(f^) 
X = Aa; 



(1) 



5ij , which are 



where Aij — &ij -\- and = dxi /dxj 
the spatial derivatives of the lensing deflection angle. $ij 
can also be written as dx^dx where <E> is sometimes called 
the lensing potential. Matrix A can be alternatively written 



Let us also define the ellipticities of the image as: 

Qii — Q22 



£2 



Qii + Q22 
2Qi2 

Qii + Q22 



(4) 



In the absence of the PSF, the quantities e\ and £2 are 
often thought to be good estimators for 71 and 72 up to the 
first order in the shear. Let us find out if they are indeed 
unbiased shear estimators. The observed quadrupole can be 
rewritten from using eq.(IT} in eq.(l2]): 



Qij = j d^xXiXjfs{A ^x) 

= !det(A)| / d2f(Af),(Af)j/s(f) 



(5) 



Note that the last step of the above equation is achieved by 
redefining A~ Keeping up to first order in k, 71 , and 

72 in eq.©, we get: 

Q11-Q22 = (l + 4ft)(0fi-0f2) + 27i(Qn+Qf2) (6) 

O12 = (l + 4fc)Qf2-f 72(Qfl+Qf2) 

O11 + Q22 = (l + 4A^)(Qfi + Qf2) + 27i(Qfi-Qf2) 

+ 472Qf2 

where Qfj are the quadrupole moments of the original galaxy 
image defined as: 



d^xXiXjfs{x) 



(7) 



Note that the two light centers defined in the image and 
source planes coincide. Based on eq.([6}, we find: 



ei = 



£2 



where 



s 



s 

£2 



4 + 271 



1 -f 271 ef + 272ef 
1 



(8) 



4 + 271 



S\2 



4 + 272 



1 + 271 ef + 272ef 
1 



4 + 272 



S\2 
£2 



r, S S 

27261 £2 



o S S 
■ 2716162 



?22 



111 + Q22 

2Qf2 
111 + Q22 



(9) 



Given that the surface brightness distribution of the original 
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galaxy image is statistically isotropic, we have (£1^2) = and 
(ef ) = 0. Therefore, we find 



(ei) = 271 
(£2) = 272 



(10) 



This result, EQ. pOf) . clearly shows that ei and £2 are not 

unbiased shear estimators, as ((ff) ) and ((£2) ) in the 
multiplicative factors depend on the galaxy morphology dis- 
tribution, and cannot be reduced to constant factors. (Also 
see Eq.(3.29) of Bernstein & Jarvis (2002), or Eq.(9.5.26) of 
Weinberg (2008).) 

One can construct an unbiased estimator of the shear, 
if one keeps three quantities from each lensed galaxy im- 
age: Qii — Q22, 2Qi2, and Qn -I- Q22, and use the ratios 
of their averages. Assuming statistical isotropy of intrinsic 
galaxy shapes in eq.©, and keeping up to first order in 
shear /convergence, we have (also see Eq.(9.5.30) of Wein- 
berg (2008)): 



1 (Qii - Q22} 

2 <Qii + Q22) 

(Q12) 
(Qii + Q22) 



71 
72 



(11) 



This form of shear estimators is not conventional, as one has 
to keep more than one quantities from each galaxy image for 
each shear component. It is this class of estimators we shall 
discuss in this paper in detail. 

One may wonder whether unbiased shear estimators in 
the conventional form ever exist. The answer is yes, at least 
when the PSF is absent. For example, we find the following 
unbiased shear estimators: 



l + £l 



4\ Vi-£i 
1 
4 



(1 



71 
72 



(12) 



Ea.(|12|l can be checked by applying Taylor expansion 
of ln[(l -I- £i)/(l - £0] (* = 1.2) to the first order in 
shear/convergence using eq.(|8}. Eq. (|12|) defines a special 
type of conventional shear estimators that are accidentally 
found by us. It is now immediately interesting to ask if there 
exist other types of unbiased shear estimators in the conven- 
tional form. We study this issue specifically in the next two 
sections. If the readers wish to go directly to the relevant 
sections on the new estimator, read on from ij3l 

For notational convenience, we shall abbreviate "con- 
ventional shear estimator" as "CSE" in the rest of the pa- 
per. Once again, by CSE we mean the shear estimators that 
are made of just one number measured from a galaxy image 
for each shear component. 



2.2.1 The Spin of the CSE 

To study the forms of the CSEs, it is useful to first consider 
their properties under coordinate rotations. Suppose we ro- 
tate the coordinates (x, y) clockwise by an angle 9. The new 
coordinates (s*,?/*) are related to the old one via the fol- 
lowing relation: 



— X cos 6 — y sin 
= X sin 9 + y cos 9 



(13) 



If we write the position vector as a complex number of the 
form X -f iy, where i is the complex unit, the coordinate 
transformation under rotation can then be written as: 



-I- iy* = (a: + \y) exp(i&) 



(14) 



For notational brevity, we shall generally use X to denote 
the value of any quantity X in the new coordinates that are 
rotated clockwise by an angle 9 with respect to the original 
coordinates. 

Let us now discuss how cosmic shears and their CSEs 
transform under coordinate rotation. The definitions of the 
shear components (shown in the beginning of ^2.1\ involve 
spatial derivatives; thus, their transformation rules under 
coordinate rotation can be found from the chain rule: 

dx" d ^dy" d 



dx dy^ 
dy'' d 
dy dy^ 



cos 9 



d_ 

dx dx dx^ 
d_ _ dx" d 
dy dy dx^ 

From ea.(|13p. we get 

d^^d^l 

dx dy 

^ = -— =sine 
dx dy 

Therefore, we have: 

d . d ( d . d\ 

d^ + 'W^[d^^%)'"''^''^ 
Taking the square of ea. (|17l) . we find: 



(15) 



(16) 



d{x0)'2 a(y»)2 

d^_^ 
dx'^ dy^ 



+ 12 



+ i 2 



dx^dy 

d' y 



(17) 



(18) 



' dxdy j 



exp{\29) 



Therefore, the shear components, which are 2nd order 
derivatives of the lensing potential $, also transform under 
coordinate rotation as: 



7i + 172 



(71 + 172) exp(i26') 



(19) 



Because of this property, we usually say that cosmic shears 
form a spin-2 quantity. In general, a complex quantity, say 
n, is called a spin-n quantity if it transforms as D* = 
nexp(inS) under a clockwise coordinate rotation of angle 



2.2 CSE in the Absence of the PSF 

In preparation for our main theme of this section, we discuss 
the spin properties of cosmic shears and their estimators in 
i]2.2.1l We then study the forms of CSEs in the absence of 
the PSF in ^2^2] 



Now let us discuss shear estimators. It is straightfor- 
ward to see that the CSEs defined in eq. (|12|l do not form a 
spin-2 quantity. More generally, assuming that Fi and r2 are 
the CSEs for 71 and 72, respectively, then, unlike 71 -I- 172, 
Fi -f ir2 is not necessarily a spin-2 quantity. However, it 
turns out that we can regularize any CSEs by turning them 
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into components of a spin-2 quantity. We have the following 
lemma: 

Lemma 2.1. Based on any pair of CSEs (ri,r2), one can 
build a new pair of CSEs (r'i,r2) to form a spin-2 quantity 
through the following procedure: 



1 Z"^" 

r'l + iFa = — / dd expi-m) (r? + ir: 

^'^ Jo 



(20) 



Proof. Firstly, under a clockwise coordinate rotation by an- 
gle 9o, we have: 



rf " + irj"" 

1- J\eeM~m (r?+'' +ir^°+*) 



2^ 



(21) 



d6l'exp[-i2(6'' -^o)] (rf 



ir? 



= (r'l + ir^) exp(i2eo) 

Therefore, r'l + ir2 form a spin-2 quantity. To show that 
(ri, r2) are a pair of CSEs, let us take the ensemble average 
on both sides of eq. (|21l) : 



(r'l+ir^) 

d6lexp(--i26')(r? + ir2) 



(22) 



1 

2^ 



Since by definition, Ti and r2 measure the shear values in 
the rotated coordinates, we have. 



(23) 



^r?^ = 7i = 7i cos 26 ~ 72 sin 28 

{^2} = 72 = 71 sin 261 + 72 cos 29 
In a more compact form, we can write ea. (|23p as: 
(r? + ir^) = (71 + 172) exp(i26l) (24) 
Using eq. ((24l) in eq. we get: 

(r'l+ir^) =7i+i72 (25) 

which proves that (r'i,r2) are indeed also a pair of CSEs. 
□ 

Due to the invariance under a coordinate rotation of 
angle 27r, one can always decompose any shear estimator 
into components of integer spins. Ea. (|20|l essentially defines 
a procedure of isolating the spin-2 components of any CSEs 
using Fourier transformation. Since the cosmic shears form a 
spin-2 quantity, only the spin-2 components of any CSEs are 
the relevant/principle components of the estimators. This 
point is further supported by the fact that the ensemble 
averages of spin-n (n 7^ 2) components of any CSEs are 



zero, i.e.. 



1 Z"^" 

— J rf6lexp(-ine)(r?+ir^) 



(26) 



27r 



d6 exTp{—m9) (71 -I- 172) exp(i2^) 



= ( if n / 2 ) 

Therefore, we only need to focus on spin-2 CSEs from here 
on. 



2.2.2 Spin-2 CSEs 

In the weak lensing limit, i.e., when the cosmic shear pa- 
rameters (71,72,^) are small, any spin-2 CSEs (ri,r2) can 
be Taylor expanded to the first order in shear as follows: 

ri(7i,72,'^) (27) 



(ri)o + 7i +72 



d-fi 



dj2 



\ dn J 



r2(7i,72,K) 

- '-..-(Sf).-'(£).-(t.. 

where {X)o {X is any quantity) denotes the value of X at 
71 = 72 = K = 0. Since (ri)o, (r2)o, {d^Ti)o, and (9Kr2)o 
are all spin-2 quantities, their ensemble average must van- 
ish. On the other hand, the coefficients associated with 71 
and 72 in ea. (|27|l can be decomposed into spin-0 and spin-4 
components as follows: 

(ri(7i,72,K)) =7i(^ + Bi) +72(C4-B2) (28) 
{r2(7i,72,'i)) =71 (B2-C)+72 (A-Bi) 



where 



A : 



(a^iri + a^,r2) 



(29) 



- 972r2) 



B2 = -(9^,ri + 9^,r2) 

As shown in Appendix A, A is a scalar, C is a pseudo scalar, 
B\ -|- i_B2 is a spin-4 quantity. The ensemble averages of B\ 
and B2 must vanish. The ensemble average of C vanishes 
if galaxy images have parity symmetry along any direction 
in the plane of the sky statistically, which is assumed to be 
true in this paper. Consequently, for (ri,r2) to be spin-2 
CSEs, we only require (A) = 1. This actually implies that 
for any individual galaxy, A — 1. The reason is that A of 
any single galaxy does not change under coordinate rotation 
of random angles, and is equal to (A) because the galaxies 
generated by rotations of a single galaxy form a complete 
set of statistically isotropic samples (i.e., there are no spe- 
cial directions). As a result, any spin-2 CSEs (ri,r2) must 
satisfy the following necessary condition: 



dVi ^ dr2 

dfi 972 



(30) 
1 directly 



This is also a sufficient condition, because A 
implies that (A) — 1. 

In general, the CSEs are functions of a certain num- 
ber of shape parameters {e.g., the multipole moments of an 
image). The functions can be very complicated, and are cer- 
tainly not fixed by the requirement given by ea. (|30p . How- 
ever, for galaxies whose shapes are described by only three 
parameters (perfect ellipses), any CSE should reduce to a 
function of just three variables. In this case, we find that 
any pair of spin-2 CSEs must reduce to a unique form, which 
is sufficient for us to judge whether CSEs are convenient in 
practice: namely, if we find that the resulting form is highly 
non-linear even for such a simple case, then it is reason- 
able to conclude that CSEs are not so useful for accurate 
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shear measurements from more reahstic galaxy shapes as 
well. This is shown in the rest of this section. For clarity, 
we refer the readers to Appendix B for the mathematical 
details/proofs for some of the statements made hereafter in 
this section. 

Let us consider a set of galaxies whose surface bright- 
ness profiles can be parametrized as fs{R) with R = a{x^ + 
y^) + h{x^ — y^) + 2cxy, where x and y are the coordinates, 
fs{R) is a function of a fixed form, and (a, b, c) are the three 
parameters determining galaxy shapes. For the images to be 
ellipses, we require the following three things: 1. fs{R) de- 
cays sufficiently fast when R becomes large; 2. a -I- fe > 0; 3. 

- 6^ > c^. For example, if fs{R) = H{Rc - R) {H is the 
step function) and (a, b, c) satisfy the above conditions, the 
galaxy surface brightness is then distributed evenly inside 
the ellipse defined by a(x'^ + y^) -\- b{x^ — y^) + 2cxy ^ Re- 
When such images are weakly lensed, the three conditions 
are not violated, and {a,b,c) becomes {a',b',c') without 
changing the form of fs- In other words, weak lensing does 
not introduce additional degrees of freedom to the galaxy 
shapes. Note that otherwise, one has to consider using more 
than 3 parameters to construct shear estimators. Among 
the three parameters, there are indeed only two degrees of 
freedom useful for shear measurement: the ratios of the pa- 
rameters. This is because the overall amplitudes of (a, b, c) 
only change the galaxy size, not its shape. As shown in Ap- 
pendix B, the ellipticities (ei, £2) defined in H2.1l are directly 
equal to {—b/a, —c/a), therefore, we can write the shear es- 
timators as functions of only ei and €2- 

We can further show that (Fi, F2) must take the follow- 
ing form: 

Fi +iF2 = (ei +ie2),gH (31) 

where u = €\ + e\, and p is a one- variable complex function, 
whose form is to be determined later in this section. To 
see why eq. (|3ip is true, one can use the Taylor expansion to 
write Fi and r2 as power series of ei and £2 Q. Consequently, 
one can write Fi -|-iF2 as power series of ei -|-ie2 and ei — ie2, 
whose spins are 2 and —2 respectively. Since Fi + iF2 is 
a spin-2 quantity, each term in the power series must also 
be a spin-2 quantity. Therefore, in each term of the power 
series, the power on ei -f ie2 must be larger than that on 
ei — ie2 by exactly one, i.e., each term must take the form 
of a(ei -f ie2)(ei + ^2)" ) where n is any non-negative integer, 
and Q is a coefficient which can be any complex number at 
this point. As a result, the shear estimators must have the 
form defined in eg. pijl . To find out what g{u) is, let us write 
it as giiu) -\- ig2{u). Eq. pi|l then becomes: 

Ti = tigi{u) - e2g2{u) (32) 
F2 = ei5'2(M) + £251 (w) 
Using the constraint in eg. pop . we find: 

2 ^ f^i + l^ (33) 
071 C72 

m\de2_ dl\de2_ dr2dej_ ^Fa de2 

dei 971 de2 d^i dei 8^2 dt2 972 

^ We do not consider shear estimators that cannot be Taylor 
expanded as power series of the galaxy shape parameters in this 
paper. Such shear estimators likely suffer numerical instabilities 
in practice. 



= 2(2 - u)gAu) + Au(l ~ u)'^ 

du 



It is interesting to note that ea. (|33p does not place any con- 
straints on g2{u), i.e., it can be any real function. This is 
because 52(1*) simply adds unnecessary even-parity terms 
into odd parity ones, and vice versa, without affecting the 
ensemble averages and the spin of the shear estimators. For 
convenience, we set (72(1*) ~ hereafter. 

Ea. (|33|l is a typical first-order ordinary differential 
equation. It can be solved by introducing an integrating fac- 
tor k{u) which satisfies: 



k{u){2 -u) = 4- [2k(u)u(l - u)] 
du 

Multiplying both sides of eQ. (|33|) with k{u), we get: 

4- [2k(u)un - u)gi{u)] = k(u) 
du 

It is now straightforward to solve both ea. p4|) and eg. ((35 
The results are: 



(34) 



(35) 



k{u) oc (l-ii)~^/^ 



(36) 



giiu) = - (l + CVT 



where C is a real number constant. To guarantee that Fi and 
F2 do not diverge when ei and £2 approach zero, we need 
C = —1. Finally, we find the unique form for the spin-2 
CSEs: 



Fi + iFa = (ei -|- ie2 



f2 4_ f2 
^1 + ^2 



(37) 



Regarding the uniqueness, it is useful to note that if we 
transform the CSEs defined in eo. (|12|) into spin-2 shear esti- 
mators using the procedure given in eg.(|20[l. we achieve the 
same shear estimators as those shown in eo. (|37|) . 

We have shown that the principle components (spin-2) 
of any pair of CSEs have to take specific and highly nonlinear 
forms for galaxies of elliptical shapes. This feature makes 
CSEs not convenient in practice (e.g., in the presence of 
noise) . 



2.3 CSEs in the Presence of the PSF 

Any CSEs which correct for the PSF effect also have to 
reduce to the forms given in eg. pzp in the limit of zero PSF 
size when the galaxy images have pure elliptical shapes. For 
this reason, the conclusion in the previous section is already 
sufficient to argue against the usefulness of CSEs in practice. 
For academic interests, we provide the following arguments 
for why CSEs may not even exist in the presence of the PSF: 
In the presence of the point spread function, structural 
details of galaxy images on scales smaller than the size of 
the PSF are smeared out. This implies that there are only a 
finite number of shape parameters (e.g., multipole moments 
up to some order) available for constructing shear estima- 
tors. On the other hand, the derivatives of the lower or- 
der shape parameters (e.g., lower order multipole moments) 
with respect to the cosmic shears depend on the higher or- 
der shape parameters because of the PSF, suggesting the 
requirement for an infinite number of shape parameters to 
form the shear estimators. Combining the above two reasons, 
we find it unlikely to form CSEs when a PSF is present. The 
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mathematical details of the above statements are given in 
Appendix C. 



{difod2.fo) 



{(drfor + {d2foV + A) 



= -72 



where 



3 A NEW WAY OF ESTIMATING SHEARS 

As searching for optimal shear estim a tors is active l y on- 
going nowadays (iHevmans et all I2OO6I : iMassev et aH I2OO7I : 
iBridle et al.ll2009l . 120101 ). it is important to realize that CSEs 
("conventional shear estimators," by which we mean the 
shear estimators that are made of just one number from a 
galaxy image for each shear component) are hard to use in 
practice due to their unavoidable complex forms even in the 
absence of the PSF (simpler forms, such as the quadrupole 
moments, are biased estimators, as shown in § 2.1). 

Therefore, existing shear estimators of the conventional 
type must quantify the bias factor when estimating the 
shear, which can be achieved numerically (se e, e.g., Erben et 
al. 2001, Bacon et al. 2001 , or most recently. i Hevmans et all 
I2OO6I : iMassev et all l2007l : iBridle et all bOld ) or estimated 
analytically {e.g., shear susceptibility in KSB [Kaiser et al. 
1995] and derived methods, or responsivity factor in Bern- 
stein & Jarvis 2002 and similar methods), although most 
people have been mainly focusing on the systematic er- 
rors caused by the photon noise and the PSF. However, to 
achieve percent or even sub-percent level accuracy in cosmic 
shear measurements, it does not seem enough to completely 
rely on numerical tests using computer-generated galaxies 
of limited morphology richness, or approximate analytical 
methods. Unfortunately, in the presence of PSF, most of the 
existing shear measur ement methods are too complicated or 
too model-dependent (|Voigt fc Bridlell2010l : iBernsteinllioiol ') 
to allow for an accurate analytic analysis of the systematic 
errors in their shear estimators. 

The method of Z08 (see also lZhandbOlOal for the treat- 
ment of photon noise and the pixelation effect) is easily 
amenable to the corrections described in ea. (|10|l . and can 
also account for the PSF correction. Not only is it simple, 
but also well defined regardless of the morphologies of galax- 
ies and the PSF. We show here how to properly use this 
method (instead of using it as CSEs) to recover the cosmic 
shear in an unbiased way. 



3.1 The Idea 

The basic idea of Z08 is to use the spatial derivatives of 
the galaxy surface brightness field to measure the cosmic 
shears. It relies on the fact that gravitational lensing does 
not only distort the overall shape of the object, but also 
locally modifies the anisotropy of the gradient field of the 
surface brightness. As it allows for using the shape infor- 
mation from galaxy substructures, the method of Z08 can 
potentially improve on the signal-to-noise ratio of the shear 
measurements. 

It is shown in Z08 that the shear measurement should 
be carried out in the Fourier space, in which any PSF can 
be transformed into the desired isotropic Gaussian form 
through multiplications, and the spatial derivatives of the 
surface brightness field can be easily measured. The cosmic 
shear can be estimated using the following relations: 



A = ^V/o • V(VVo) 



(39) 



1 {{d^fof - {d2fof} 

2{{d^for + {d2for+A) 



P is the scale radius of the isotropic Gaussian PSF W0, which 
is defined as: 



■ exp 



"2/32 



(40) 



fo is the surface brightness field, di denotes d/dxi. As shown 
in Appendix D, the method of Z08 effectively utilizes the 
quadrupole moments in the Fourier space to measure the 
cosmic shears. 



3.2 A New Unbiased Estimator 

Now, here is an important point: in order to implement this 
method, we must make it clear what we mean by the angular 
brackets in eq.((38|. First, we need to measure the derivatives 
of the surface brightness and average them within a single 
galaxy. Let us denote this averaging by ()g, and write: 

1 {{difof - {d2fof), 
2{{difoY + [d2foY + A), 

{difod2fo)g 

{{d,foY + {d2foY + A), 

Of course, these are still extremely noisy as they use only 
one galaxy. The question is then, "how do we average these 
quantities over many galaxies to obtain an unbiased estima- 
tor of the shears?" 

If one uses these quantities as if they were the CSEs, 
then one would simply average them over many galaxies. 
However, this will produce a biased estimator: 



1 {{d,fof - {d2fof), 

2{{d^foY + {d2foY + A), 
{difod2fo)a 

{{d,foY + {d2foY + A), 



-71(1 -<5i) (41) 



= —72(1 — &) 



where 5\ and 82 are the ensemble averages of functions of 
multipole moments of the galaxy images in Fourier space, 
and {)en. denotes the ensemble average over many galaxies. 
The derivation of the forms of 5i and 82 is given in Appendix 
D. (5i and 82 are generally nonzero and dependent on the 
galaxy morphology. 

Instead, we need to take the ensemble averages of the 
numerator and the denominator separately first, and then 
divide them to obtain an unbiased estimator: 



1 {{{dlfo? - {d2fof},)en 

2{{{difor + {d2foy + A)g},„ 

{{dlfod2fo)g)en 



= -71 



(42) 



= —72 



-71 



(38) 



{{{difor + {d2for+A)g},,, 

This is the main result of this paper, and the form of the 

2 This bias was not detected in the numerical calculations pre- 
sented in Z08, as the number of galaxies (10'*) used for the average 
was still too small. 
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unbiased estimator that we propose to use for the actual 
analysis of the weak lensing data. 

Of course, one could divide the left hand sides of ea. (|4ip 
by 1 — <5i and 1 — 52 to obtain an unbiased estimator. This is 
similar to correcting the measured shear for a multiplicative 
bias that is evaluated from the same ensemble of galaxies. 
In this sense, eq. (|4ip provides the exact definitions for the 
multiplicative biases for 71 and 72. However, since 5i and 
S2 in eQ. (|41[l involve many high order Fourier-space multi- 
pole moments of the surface brightness field, evaluation of 
these terms from simulations (which are incomplete anyway) 
can be highly uncertain. Even worse, the multiplicative bias 
mentioned here is not even a constant, but depends on the 
morphological distribution of the galaxies. This makes the 
conventional way of measuring shear correlation functions 
even more challenging, as one must take into account the 
correlations of the multiplicative biases, as will be shown in 

m 

In summary, according to eg. 1)42 1) . for each shear compo- 
nent, two quantities from each galaxy should be kept, and 
the ratios of their ensemble averages yield unbiased esti- 
mates for the corresponding shear components. Finally, it is 
important to note that, to efficiently use eg. (| 42 [I . the sur- 
face brightness of each participating galaxy should be nor- 
malized to have roughly the same maximum value, so that 
faint galaxies are not much less weighted than their brighter 
counterparts. The details regarding the optimal weighting 
scheme as a function of the galaxy luminosity should also 
take into account the photon noise. This is a separate topic, 
and will be studied in a future work. 



3.3 Comments on Errors due to Finite Number of 
Galaxies 



Strictly speaking, Eq.((42|) holds when we average over an 
infinite number of galaxies; however, as we shall show in this 
section, the error that we make by having a finite number 
of galaxies for averaging is much smaller than the statistical 
errors, and thus the estimator remains unbiased for practical 
applications. 

For simplicity, we use eg. (|ll[) rather than our main 
equation [eQ. (|42p ] in the following discussion, but the con- 
clusion will be the same for eg. (|42p . 

Let us use (X) n to denote the average of the quantity 
X over galaxies. From eq.((6|, we get: 



(Qll-Q22)^ = (l+4Ac)(Qfi-Qf2)^ 

+ 27i(Qfi + Qf2)^ 

(Ql2)^ = (l+4^.)(Qf2)^+72(0fl+0f2) 

{gii + Q22)^ = (l + 4^.)(Qfi + Qf2)^ 

-f 271 (Qfi - Qfa)^ + 472 (Qf2)^ 
Consequently, we have: 

1 (Qll — Q22)ff 

2 (Qii -f Q22)^r 

(<3l2)jv 



(43) 



^Ai +71(1 - A?) -72A1A2 (44) 



{Qll+Q22}j 



-72(1 



- 71A1A2 



where 



Ai = 



A, = 



(Qfi-Qf2)^ 
{Qii + Q22) jv 

JV 



(45) 



Here, the terms Ai and A1A2 contribute to random errors 
because their ensemble averages vanish, whereas the terms 
Ai and A| lead to systematic biases because their ensemble 
averages do not vanish. Fortunately, as Ai,2 scales as l/\/iV 
and 71,2 ^ 1, the amplitudes of such systematic biases are 
always much smaller than the sizes of the statistical errors. 
Therefore, the results from this new type of shear estima- 
tors may be regarded as unbiased for practical applications. 
Numerical verifications will be given in 3S1 



SHEAR STATISTICS 
CORRELATIONS 



N-POINT 



The cosmic shear field can only be probed statistically. This 
is mainly due to the intrinsic variations of the galaxy shapes 
and the spatial fiuctuations of the shear components. As a 
result, the shear statistics is usually studied in the form of 
n-point spatial correlation functions of the shear field. The 
previous discussions and measurements in the literature are 
based on "conventional" shear estimators (CSEs), i.e., one 
often assumes that the following is true: 



{r)e„ = 7 



(46) 



where 7 can be either 71 or 72, and P is a CSE for 7. For 
individual galaxies, eg. (|46p implies: 

r = 7 + * (47) 

where 'if satisfies {^')en ~ c|f|. It is usually assumed that 
^f's of different galaxies do not correlate with each otheiQ. 
Therefore, the n-point correlation functions of the shear field 
can be directly measured by the correlations of F's. 

However, in ^2.2\ we have shown that such a F at least 
does not exist in a convenient form. Instead, as proposed in 
H3.2I we can use the new form of shear estimators defined in 
eg. (|42p to probe the cosmic shear in an unbiased way. Let us 
now find out how to measure the n-point shear correlation 
functions with the new form of shear estimators. Numerical 
examples are given in ^ 

For notational convenience, the type of shear measure- 
ment in eg. (|42p can be symbolized as follows: 



: 7 



(48) 



where 7 can be either 71 or 72 , and A and B are properties of 
a galaxy, such as those defined in eg. (|42p . Similar to eg. (|47p . 
eg. (|48p implies the following: 



(49) 



^ Note tliat for our purpose, it is not necessary to know the form 
of 

* This is at least true if the relevant galaxies are separated by 
a large physical distance. Detailed discussions about the correla- 
tions of ^'s belong to the topic of "Galaxy Intrinsic Alignment", 
which is beyond the scope of this paper. 
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where C satisfies {C}en = 0. If we assume that the C of 
any galaxy does not correlate with the _B's and C's of other 
galaxie^, the n-point correlation functions of the shear field 
can be probed using the following relation: 



(7(^1)7(52). •.7(^n))^„ 



{A{xi)Aix2)...A{x„))^ 
{B{xi)B{x2)...B{x,,))^ 



(50) 



The ensemble averages are taken over a large number of 
galaxies whose relative positions Xi — Xj {i, j = 1,2, . . . ,n) 
are fixed. In practice, the n-point shear correlation functions 
can be measured using: 



(51) 



Y:,A{xi,)A{xi^)...A{x,J 
Y,B{x,,)B{x,,)...B{xiJ 

where the sum is taken over all the galaxy groups that sat- 
isfy the positional constraints. Note that the ratio is taken 
after the summations. The standard deviation (a) of the cor- 
relation function in such a measurement can be calculated 
as follows: 



J2B{xi,)B(x^^)...B{x,„ 

/ j:aHs.^)a^x.,)...aHx,j 

\[J:B{x,,)B{x,,)...B(x,jY 
N{A\x,)A'{x2)...A\x„))^^ 



(52) 



m {B{xi)B{x2) 

1 f {A^)en^ 



■B{Xn)) 



where A'' is the total number of galaxy groups {e.g., the 
number of galaxy pairs for 2-point correlations) used. 

To summarize, in the new type of shear measurement, 
the shear correlation function should be measured using the 
ratio of two ensemble averages, as shown in ea. (|50p . If B in 
ea. (|48p is viewed as a multiplicative bias, we need to measure 
the correlations of these multiplicative biases as well in order 
to get the correct shear correlation functions. 



5 NUMERICAL TESTS 

In this section, we show how accurately one can recover the 
cosmic shears and their 2-point correlation functions with 
the method proposed in SJ3] and SJJ] For a comparison, we 
also show the results using the method of Z08 in the "con- 
ventional" (but wrong) way, i.e., ea. (|41| . but without tak- 
ing into account the biases. Si and 621] Since our focus is 
to demonstrate principles, we do not include photon noise 
or the pixelation effect in this paper, and we approximate 
the PSF as an isotropic Gaussian. Including these effects is 
straightforward. (Note also that the conventional estimators 
yield biased results even in this idealized case.) 

More comprehensive tests of the method of Z08 have 



^ Such correlations again belong to the topic of "Galaxy Intrinsic 
Alignment," which is not considered in this paper. 
^ Note that Z08 did not detect these biases because it did not 
use enough galaxies (10'*) in the numerical tests. Here, we use 10^ 
galaxies (or pairs of galaxies for the 2-point correlation function) 
for each test. 



been presented in IZhand (|2010bl ). which further improves 
the accuracy of shear measurement by including the second 
order shear/convergence terms. As shown in that paper, the 
second order correction is proportional to the convergence 
K. We simply set k — in the numerical simulations here. 
The residual systematic error on the recovered shear (71 
or 72) therefore should be expected to have a magnitude 
comparable to the third order terms in shear {e.g., 7?, 7172, 
7172, 72) 

all at the la confidence level. 



5.1 Image Generation 

The mock galaxy images we use in our numerical tests are 
generated by the algorithms introduced in Z08 and Zhang 
(2010a), i.e., each galaxy is generated as a collection of point 
sources. The reason is simple: one can accurately and easily 
mimic the lensing effect by displacing the points. It also al- 
lows us to generate galaxies of complex morphologies. There 
are two types of galaxies we use in this paper: 1. randomly 
oriented regular galaxies, each of which contains an expo- 
nential disk in the galactic plane (no bulge); 2. irregular 
galaxies being made of points generated by the trajectories 
of 2D random walks. For simplicity, the PSF is always an 
isotropic Gaussian function, whose scale radius is four times 
the grid size to avoid the pixelation problem. All the lengths 
in our simulations are in units of the grid size in the rest of 
this section. The dimension of the grid is 64 x 64. 



5.2 1-Point Statistics 

As our first example, we study how accurately a single input 
cosmic shear can be recovered by a large number of mock 
galaxies, i.e., the 1-point statistics. We use the regular type 
mock galaxies as introduced in ^5.1\ Each disk galaxy is com- 
posed of ten point sources which are randomly distributed 
within a radius of 7. The intensity of a point is an exponen- 
tially decaying function of its distance to the center of the 
disk with a decay length equal to 7. The galactic disk is then 
projected onto the source plane in a random direction. For 
each input shear value, we use 10^ mock galaxies to recover 
the shear. 

To quantify the accuracy of shear recovery, we adopt 
the standard technique in the weak lensing community by 
using the "multiplicative bias" rm and the "additive bias" 
Ci, which are defined as: 



measured /i , \ input 

71 = (1 + "^i)7i + ci 

measured / -i , \ input , 

72 = (1 + "^2)72 + C2 



(53) 



Our simulations use six sets of input shear values (71, 
72). They are: (0.05, -0.05), (0.03, -0.03), (0.01, -0.01), 
(-0.01, 0.01), (-0.03, 0.03), (-0.05, 0.05). The recovered 
shear values as well as the linear fitting results for 71 and 72 
are shown in Table [1] and Table [2] respectively. Note that we 
also list th e values of y^'s an d Q's for the goodness of linear 
fitting (see lPress et al.l|l992l for details). For a comparison, 
we show in the last column of each table the quality of the 
shear recoveries using the "conventional" (but wrong) way, 
given by eQ. (|41l) . 

The tables show that the shear recovery can be very 
accurate if we use the method of Z08 in the proper way. 



© 2010 RAS, MNRAS 000.[TlfT3l 



Cosmic Shears Should Not Be Measured In Conventional Ways 9 



Input 7i 


7i measured in 
the proper way 


7i measured in 
the "conventional" way 


0.05 


U . UOUUI^U. uuuuo 


n 0^917-1-0 nnnoQ 

u . Utjz ± / ^u. uuuuy 


u.uo 


U . UOUU4I1IU. uuuuo 


u.uo -Lo^mu.uuuuy 


0.01 


0.01006±0. 00008 


0.01049±0. 00009 


-0.01 


-0.01002±0. 00008 


-0.01046±0.00009 


-0.03 


-0.02997±0. 00008 


-0.03126±0.00009 


u.uo 


— u.uouuomu.uuuuo 


— u.uoz j.oiiiu.uuuuy 


Linear Fitting 
Results 


Using 71 's from 
the proper way 


Using 71 's from 
the "conventional" way 


mi 


(0.7 ±1.0) X 10-=* 


(43.4 ± 1.1) X 10-=^ 


ci 


(1.6 ±3.4) X 10-5 


(1.2 ± 3.7) X 10-5 


ixl Qi) 


(0.58, 0.96) 


(0.58, 0.96) 



Table 1. In the middle and right columns of the upper part of 
the table, we list the measured 71 's from the proper way and 
the "conventional" way of using the method of Z08, respectively. 
The lower part of the table shows the multiplicative biases, the 
additive biases, and the goodness of the linear fittings for both 
cases. The definition of the linear fitting here is given in cq. ll53l l. 



Input 72 72 measured in 72 measured in 

the proper way the "conventional" way 

0.05 0.05011±0.00008 0.05226±0.00009 

0.03 0.02996±0.00008 0.03126±0.00009 

0.01 0.00989±0.00008 0.01033±0.00009 

-0.01 -0.01010±0.00008 -0.01055±0. 00009 

-0.03 -0.03000±0.00008 -0.03126±0. 00009 

-0.05 -0.05011±0.00008 -0.05225±0. 00009 



Linear Fitting Using 72 's from Using 72 's from 

Results the proper way the "conventional" way 

m2 (1.4 ± 1.0) X 10-^ (44.2 ± 1.1) X 10-^ 

C2 (-4.1 ±3.4) X 10-5 (-3.5 ± 3.7) X 10-5 

(xi, Q2) (3.3, 0.51) (3.4, 0.50) 



Table 2. Same as table [T] except that it is for 72. 



i.e., eq. (|42|l . On the other hand, if Z08 is used in the "con- 
ventional" way, i.e., eq. (|4ip . we find nonzero multiplicative 
biases (mi, 7712) for both 71 and 72. The additive biases ci 
and C2 are consistent with being zero. These results agree 
with our conclusions in ^i3.2\ 



5.3 2-Point Correlations 

Let us now test the recovery of the 2-point shear corre- 
lations using the method of Z08 in the way that is pro- 
posed in ^ In each test, we use two large groups of galax- 
ies that are given (71, 72) and (7!^, 72), respectively. Each 
group of each test contains lO'^ galaxies. The values of 71, 
72, 7i, and 72 vary from galaxy to galaxy. They are as- 
sumed to be normally distributed with the following co- 
variance matrix: {7?) = (7^) = (7I) = {72^) = 0.04^ 
(7172) = (7172) = (7W2) = (7172) = 0, and (717;) and 
(7272) are to be specified in each test. The purpose of the 
test is to find out how accurately the 2-point correlations 
{7171) ^-nd (7272) can be recovered. For the tests presented 
here, we use the irregular type of mock galaxies, each of 
which is made of ten point sources generated by the 2D ran- 
dom walks. Each step size of the random walks is a random 
number between and 2. The radius of each galaxy is lim- 
ited to be less than 7. 

In table [3] and [H we report the results of six tests with 
six different sets of ((717O, {7272)): (0.001,0.001), (0.0006, 
0.0006), (0.0002, 0.0002), (-0.0002, -0.0002), (-0.0006, 
-0.0006), (-0.001, -0.001). As in to characterize the 
accuracy of the method, we again use the multiplicative bi- 
ases (mil, m22) and additive biases (cn, C22) that are de- 
fined as follows: 

/ f\measured i \/ f \ input , frA\ 

(7171) = (1 + mii)(7i7i) +C11 (54) 

/ / \ measured /i ■ \ / / \ input , 

(7272) = (1 + m22)(7272) + C22 

We also show the results from the "conventional" way of 
using Z08 for a comparison. 

Our results again show no systematic errors for the 
proper way of using Z08. In contrast, the "conventional" way 
tends to underestimate the amplitudes of the shear correla- 
tions (negative multiplicative bias). The signs of the mul- 
tiplicative biases in the "conventional" cases are opposite 
to those found in ^5.2\ This is because we have used two 
different types of mock galaxies. 



6 SUMMARY 

Conventionally, in the studies of weak lensing, for each shear 
component (71 or 72), one hopes to construct a single quan- 
tity from each background galaxy image, whose ensemble 
average is equal to the true value of a component of the true 
shear field. We have shown that such conventional shear es- 
timators (CSEs) do not exist in convenient forms even in 
the absence of the PSF. 

Based on the method of Zhang (2008), we have proposed 
to measure the cosmic shear in a new way: using the ratio of 
the ensemble averages of two galaxy properties to estimate 
each shear component. (Also see § 9.2 of Weinberg (2008) for 
a similar study.) We have shown that, using both analytic 
analyses and numerical examples, the new way of estimating 
cosmic shears is unbiased, and does not contain systematic 
errors to the first order in shear at least. The new type of 
shear measurement demands shear statistics such as n-point 
correlation functions to be measured in an unconventional 
way as well, but with little additional cost. 
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Input (7i7(> 


(71 7i) measured in 
the proper way 


(71 7i) measured in 
the "conventional" way 




(0.99 ± 0.03) X 10-3 


(0.71 ± 0.02) X 10"3 


6 X IQ-* 


("6 S + S") X IQ-* 


("4 5 + n 21 X 10"* 


2 X 10-* 


(1.7 ± 0.3) X 10-* 


(1.2 ± 0.2) X 10-* 


-2 X 10-"' 


(-2.1 ±0.3) X 10-* 


(-1.4 ±0.2) X 10-* 


-6 X 10-"' 


(-6.1 ±0.3) X 10-* 


(-4.5 ±0.2) X IQ-* 


— 10-'^ 


C_n QQ + n nX\ v 10-3 


f— 72 + 021 X 10"3 


Linear Fitting 
Results 


Using (7i7()'s from 
the proper way 


Using (7i7i)'s from 
the "conventional" way 


mil 


-0.001±0.018 


-0.281±0.012 


Cll 


(-0.18 ± 1.2) X 10"^ 


(-4.1 ± 7.9) X 10"® 


(X?i, On) 


(2.5, 0.64) 


(2.3, 0.68) 


Table 3. In the middle and right columns of the upper part of 
the table, we list the measured (7i7j)'s from the proper way and 
the "conventional" way of using the method of Z08, respectively. 
The lower part of the table shows the multiplicative biases, the 
additive biases, and the goodness of the linear fittings for both 
cases. The definition of linear fitting here is given in eg. II 5411. 


Input (7272> 


(7272) measured in 
the proper way 


(7272) measured in 
the "conventional" way 


10-3 


(0.99 ±0.03) X 10-3 


(0.71 ± 0.02) X 10-3 


6 X 10-"* 


(5.8 ±0.3) X 10"* 


(4.2 ±0.2) X 10"* 


2 X 10-* 


(2.0 ±0.3) X 10-* 


(1.3 ± 0.2) X 10-* 


-2 X 10-* 


(-1.7 ±0.3) X 10-* 


(-1.3 ±0.2) X 10-* 


-6 X 10-* 


(-6.0 ±0.3) X 10"* 


(-4.4 ±0.2) X 10"* 


-10-3 


(-1.03 ±0.03) X 10-3 


(-0.73 ±0.02) X 10-3 


Linear Fitting 
Results 


Using (7272)'^ from 
the proper way 


Using (7272)'s from 
the "conventional" way 


m22 


0.000±0.018 


-0.284±0.012 




(-0.61 ± 1.2) X 10-5 


(-6.2 ± 7.9) X 10-6 


ixh' O22) 


(1.8, 0.77) 


(1.6, 0.82) 


Table 4. Same as 


table [S] except that it 


is for (7272)- 
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APPENDIX A 
SPIN-2 CSES 



THE DERIVATIVES OF THE 



Under a clockwise coordinate rotation of angle 9, the cosmic 
shear components transform according to the following rule: 



7i + 172 = (71 + 172) exp(i26') 
Using the chain rule, we then find: 



d . d 



— + i— 



exp(i26') 



(55) 



(56) 



On the other hand, we know that the spin-2 shear estimators 
(ri,r2) transform as: 
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r? + ir^ = (ri + ir2)exp(i26») 

Apply ea. (f56|) onto eq. (fST)) . we get: 



Cosmic Shears Should Not Be Measured In Conventional Ways 11 

(57) / a2 - 62 



+ 1 



d'yi 



d^2 



+ i 



972 971 



(58) 



exp(i4^) 



Ea. (|58p clearly shows that the derivatives (d^f^V\ —d^f^V-^) + 
i(9^2ri + d^^V^) form a spin-4 quantity. We can also apply 
ea. (|56p onto the complex conjugate of eq. (|57p . and get: 



97? ^ 97I 

STi ^ 9r2 
971 972 



+ i 



+ i 



972 



9r| 

97? 

9r2 
971 



(59) 



Therefore, d^^V\ +9^2r2 is a scalar, and 9-^2 Fi — 9^ir2 is a 
pseudo-scalar (it has odd parity). 



APPENDIX B - SOME MATHEMATICAL 
DETAILS REGARDING GALAXIES OF 
ELLIPTICAL SHAPES 

The surface brightness profiles of the galaxies used in i]2.2.2l 
are parametrized as fs{R), where R — a(x^ + y^) + blyX^ — 
y^) + 2cxy, and x and y are the coordinate variables. Curves 
of constant values of R can be ellipses, hyperbolas, parabo- 
las, or lines depending on the values of a, b, and c. For our 
purposes, we only need ellipses. This requires a -I- fe > 
and — b'^ > according to the matrix theory. Since 
linear coordinate transformations do not spoil these rela- 
tions, weakly lensed images of these galaxies are still el- 
lipses. For the same reason, one can also easily show that 
the lensed images are still parametrized by the same func- 
tion fs{R) with R = a'(x'^ +y'^) + b'{x' -y'^) + 2c'xy, where 
a' = a(l - 2/t) - 2716 - 272c, b' = b(l - 2k) - 271a, and 
c' = c(l - 2k) - 272a. 

Let us now show that the ellipticity parameters (ei, £2) 
for galaxies of the form fs \a{x'^ + y^) + b(x^ — y^) + 2cxy^ 
are equal to {—b/a,—c/a). According to the definitions in 
ij2.1l we have: 



(fxxiXjfs \a{x^ + y^) + b{x^ — y^) + 2cxy^ (60) 



This integration can be carried out using the following linear 
coordinate transformation: 



y - 

where 



Va + b 
1 



(/3a;' + ay') 
(ax + Py') 

c 



(61) 



62 



62 



/3 = 



Va2 62 
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- 62 



- 62 



(62) 



' — 62 — c2 

Using {x',y') defined in ea. (|6ip to replace {x,y) in eq. (f60 
we get: 



(63) 



Va2 - 62 

6(a^+/3^) ,^,2 
a2 -62 

Qii + Q22 



a{j3'^ - a^) 



[x +y ) 



4bal3 



ix''- 



62 



X y 



fs(x +y ) 



62 



„ 2\ 

'{x 



62 



+ 



a{a^ + P^) ,2 12-, , ^aaP , , 



62 



r / 12 , I2\ 

fs{x +y ) 



2Qi2 



62 



+ 



2(q"+/3") , , 
— T y 
^^2— 62 



2a/3 



a" — 62 
12 , 12 



12 , l2\ 

X +y ) 



r (12 I f2\ 

fs{x +y ) 



Given that the function fs only depends on x''^ + y'^, we 
have: 



,2 7/ /2 /2\ r / 12 , 12 \ 

d x'(x - y )fs{x +y ) 



d x'x'y' fs{x' -I- y' ) = 

It is now straightforward to calculate (ei,e2): 

Qii - Q22 _ & 
" a 



(64) 



£2 



/II + (^22 
2Qi2 



(65) 



Qii + 1 



APPENDIX C - WHY OSES DO NOT LIKELY 
EXIST IN THE PRESENCE OF PSF 

For technical convenience, let us work in Fourier space. Ac- 
cording to the notations in ^2.11 we use /!,(«) and fs(k^) 
to denote the Fourier transformations of the lensed galaxy 
image fL{x'^) and the original galaxy image fs{x^) respec- 
tively. Their relations are given by the following equations: 



fs{k^) 



j2^L ik^-x^ 

axe 
axe js(x ) 



(66) 



Using the relations defined in eq.((T}, we get: 



/L(fc^) 



,2 

d X 



ldet(A)l 



det 



fdx^ 
[dxS 



(67) 



,2^S i{Ak^)-x^ 

a X e^ ' 



|det(A)|/s(Afc^) 
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Eq. (|67|) simply means that under lensing, the Fourier trans- 
formation of the galaxy image is changed from fs{k) to 
|det(A)|/s(Afc). 

Due to the presence of the PSF, the Fourier transfor- 
mation of the observed image fo is related to that of the 
lensed image /l via: 



fo{k) = W{k)fL{k) 



(68) 



where W{k) is the Fourier transformation of the PSF. With- 
out loss of generality, in the rest of our discussion, we 
use the isotropic Gaussian PSF, i.e., W{k) = VF^(fc) = 
exp(— |fc| /2). The advantage of working in Fourier space 
is that the PSF is included as a multiplicative factor, rather 
than a convolution as in real space. Combining ea. (|67p and 
ea. ([68l) . we get: 

~ ~ - (69) 



/o(fc) = W^,3(fc)|det(A)i/s(Afc) 



Since the PSF profile in Fourier space typically falls off 
quickly when the wave number exceeds the inverse of the size 
of the PSF, it strongly suppresses the power of the observed 
images on small scales. Therefore, only a finite number of 
Fourier modes are available for providing shape information. 
Recovering information on arbitrarily small scales is never 
feasible in practice due to noise and numerical problems. 

To form CSEs, let us use the multipole moments of 
galaxy images to represent the shape information, which are 
defined as: 



(fkk\k{fo{k) 



(70) 



where i and j are non-negative integers. Note that due to the 
finite degrees of freedom of the shape information, one can 
equivalently choose other basis {e.g., shapelets) to study the 
same issue without affecting the conclusion. For simplicity 
but without loss of generality, let us only consider galaxies 
that are invariant under the parity transformation x — >■ —x. 
Note that weak lensing does not change this property. For 
this type of galaxies, the imaginary part of foik) is always 
zero, and Mij's axe real. Furthermore, Mij is zero when i + j 
is an odd number. In this case, the shear estimators Fi and 
r2 can be written as functions of the Mij with i + j being 
even numbers only. For some of the lowest order Afij's, we 
can find out how they transform under lensing using eQ. (|69p : 



Mi. 



cFkkiki^Wf,{k)\det{A)\fs{Ak) 
d^k{A-^k)\{A~^ky.^Wfi{A~^k)fs{k) 



(71) 



The last step in the above equation is achieved by redefining 
Ak as k. By keeping the terms in matrix A up to the first 
order in shear, one can straightforwardly show the following: 



M20 - M02 

= ( 1 - 2«:) ( a4 - ) + «:/9' (Mfo - ivi^ 

+ 272/3' (Afg^l - Mfg) 

- -|^[2[M|a + M|2)-P^{M: 



+ M^i-2Mi2)] 



(72) 



Mil 
(1- 



2fi:)Mn 



- 72 



(mIo + Mi2 ■ 



2P^Ml2) 



+ 



Af2() + Ma 



(1 - 2«:)(a4 + a4) -f (Affo + 2*4 + A^ol) 

71 [2(Af2l, - a4) - f {Mio - Mi,)] 

72 [4Mfi - 2/3' {Mi^ + A/fs)] 



where 



Af5 = I d'kk{kiWp{k)fsik) 



(73) 



The conventional cosmic shear estimators (ri,r2) are 
functions of Mij. According to our discussion in H2.2.2I the 
two functions have to satisfy the following relation: 

1^ + 1^-2 (74) 
071 072 

In the presence of PSF, we are unable to find out whether 
one can find CSEs that satisfy the most general requirement 
given in ea. (|74|) . However, we can show that there do not 
exist CSEs satisfying a slightly stronger condition: 

^Fi dTj 

dji d"/2 

In addition to eq. H74|) . ea. (|75p simply imposes another re- 
quirement that the spin-4 parts of the derivatives of the 
shear estimators with respect to the shears are zero. We can 
rewrite ea. (|75|) for, e.g., Fi, using the chain rule as: 



= 1 



(75) 



dVi dMi, 



(76) 



Since there are only a finite number of multipole moments 
available for constructing the shear estimators, we assume 
the maximum value of i -f j is (N is an even integer), i.e., 
we have: 



1 = 



E 



dMij 971 



(77) 



Because the right side of ea. (|77p is evaluated at 71 = 72 = 
K = 0, both dVi/dMij and dMijjd^x are functions of only 
Mfj^s, I.e., the multipole moments in the absence of lensing. 

As in ea. (|72p . one can show in general that dMij/d^i 
involves higher order multipole moments, i.e., M^j,'s with 
i' + j' > i + j. In particular, when i + j — N , dMij/d^i 
depends linearly on Mfij,'s with i' + j' > N. To satisfy the 
constraint in ea. (|77|) . the coefficients in front of the terms 
proportional to Mfij, 's (i' -I- j' > N) must vanish because 
they are independent of the multipole moments with i + j^ 
N. As a, result, we find that dVi/dMij has to vanish when 
i + j = A'^. In other words, we have: 



E 



dVi dM,j 
dMij d-yi 



(78) 



We can now recursively use the above reasoning to show that 
there does not exist an A*' that can satisfy ea. (|77|l . There- 
fore, we can never find CSEs of the type defined in eq. (|75p . 
Though, one may still expect to find CSEs based on the 
most general requirement defined in ea. (|74|) . However, even 
in this case, CSEs must reduce to highly nonlinear forms 
for galaxies of pure elliptical shapes when the PSF effect is 
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small, as we have shown in i)2.2.2l This feature is already 
sufficient for arguing against their usefulness in practice. 



APPENDIX D - DERIVATION OF THE 
MULTIPLICATIVE BIASES RESULTING FROM 
A MISUSE OF Z08 



Let us now calculate the terms Si and 62 defined in eQ. (|4ip . 
The averages of the spatial derivatives of the surface bright- 
ness field of a single galaxy can be related to the Fourier 
modes of the image. The Fourier transformation has been 
defined in ea. H66p in Appendix C. Following the notations 
of Appendix C, we find: 



di.fo{x) = 
and 



(27r)2 
{difo{S)djfo{x)) 



(_ifc,)e-''=-7o(fc) 



(79) 



(80) 



-difo{x)djfo{x) 



Jgd^x J dFk J (fk' 
~ (2^0^ 
1 Jd'^kJ cfk' 



fo{k)fo{k') 



S (27r)4 
1 / (fk 



{-k,k',){2nySl{k + k')fo{k)fo{k') 



S {2ny 



fo{k) 



where S is the total area of the map containing the galaxy. 
Similarly, one can derive the following relation: 



(V/- V(vVo))^ 



1 f (fk 1^,4 
\k\ 



foik) 



(81) 



S (27r)2 

Ea. (|80p and ea. (|8ip allow us to transform eq.(j4T} into its 
version in Fourier space: 



P20 + P02 - 
Pii 

P20 + P02 - 



= -7i(l-5i) 
-72(1 — 



(82) 



where 



P.-Jd^kklki 


fo{k) 


2 


(83) 


Dn = J d^k\k\" 


ib(fc) 


2 





Note that D4 — P40 + 2P22 + ^04- It is now clear that the 
method of Z08 basically utilizes the quadrupole moments in 
the Fourier space to measure the cosmic shear. 

Using eq. H69p , we can find out how Pij transform under 
lensing: 



The last step in the above equation is achieved by redefin- 
ing Afc as k. For the isotropic Gaussian PSF, Wp{k) = 
exp(— |fc| /2). By keeping the terms in matrix A up to 
the first order in shear, one can straightforwardly show the 
following: 



P20 

nS 



P02 



(85) 



Pii) + 472/3' (Pii - PiV 



271 [pio + Pi 

Pu 



/3' iP. 



" + Pii 



2Pi2)] 



pS 

^11 



■ 72 (Pio + Pi2 



4/3' Pi2) 



+ 2kP^ (pI + Pfs) 



+ 



271/3' (pii 



Pis) 



271 [Pio - Pi2 - /3' (P4^o - Pii)] 



P20 + P02 
pio + pi2 
472 [Pi - /3' (pii + Pis)] + 2kp'd. 

Da 



d! 



2k [d! 

(Pio 



- P^Di) 



Pii) 



271 [2 
472 [2 (pl + Pi) 



(pio + pi2 



pS 
-'24 



pS 

-foe 



)] 



where 



d kk-^ 



VK8(fc)/s(fc) 



Wpik)fs{k) 



(86) 



By keeping the terms up to the first order in shear, and 
using the fact that the ensemble averages are taken over 
statistically isotropic galaxy samples, one can now directly 
find expressions for Si and S2 of ea. (|4ip : 



5i 



(pS psy 

{-t^20 ~ ^02 ) 



{Pi, + Pi2~P^Df/2Y 



(87) 



(Pio ■ 



-Pq^) (-P4^ ~ Poi) 



Pii) (pio + Pi2 



(Pio 



pS 
24 



Pii) 



{Pi, + pi2-P^Df/2Y 
(Pio + Pi2-P'Di/2y 



4(Pii) 



Pii (pji + Pis) 
{Pio + Pi2-P'D!/2f 



d kk-^ 



W,5(fc)|det(A)|/s(Afc) 



X Wf3{A-^k)fs{kf 



(84) 



+ 2/3^ 



pfi {pj + 2P3i + pg) 

{Pio + Pi2-FDi/2)' 
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